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Fig. 4 Measured probe descent module spin performance compared 
with a posteriori spin assessment. Two sets of numbers correspond to 
descent time and spin rate in the first and last data acquisition intervals. 


spin response. This assessment was derived from further model cal¬ 
culations based on drop and wind-tunnel tests. 

Five minutes past entry, the spin rate was measured as ±33 rpm 
(33.48 rpm at 4.9 min). A cursory comparison with Fig. 3 suggests 
that the spin direction at aeroshell jettison was positive, i.e., not 
reversed by ablative grooving, although insufficient data during this 
initial descent period exist to conclude this with certainty. Ancillary 
data, however, confirm that the intent of the transient spin require¬ 
ment was adequately met. Because the Orbiter-mounted relay radio 
receivers managed to lock onto the probe’s modulated signal much 
sooner than expected (~35 s instead of the specified 1-min maxi¬ 
mum), it is clear that any high-rate probe spin that might have been 
induced by entry forces was quickly damped. If this were not the 
case, the phase error of the received signal—predominantly due to 
spinning effects during initial signal acquisition—would have been 
too high to permit rapid signal acquisition. 9 

For the remainder of the descent mission, the spin assessment 
appears to provide a more conservative, lower-bound spin estimate, 
matching the flight data to within about 5 rpm. However, the as¬ 
sessment itself is based on an equation of motion, which assumed 
that the parachute was nonspinning and, hence, would not be ex¬ 
pected to accurately replicate the actual flight condition of a spinning 
parachute/spinning DSM whose composite rolling moment and roll 
damping coefficients were never measured or analyzed. Further¬ 
more, the assessment ignored the possibility of dynamic coupling 
of pitch and yaw disturbances into the probe spin. Again, this was 
because aerodynamic test and design efforts were focused toward 
meeting the transient spin requirement. For the balance of the de¬ 
scent mission (beyond 5 min past entry), it was always deemed suf¬ 
ficient that the probe simply have a nonzero spin. (It is interesting to 
note that the terminal spin rate measured during the aforementioned 
test in Earth’s atmosphere was within 4% of the final Jovian spin 
measurement.) The last point recorded, at 50.77 min after entry, 
corresponded to a spin of 14.2 rpm. 

Conclusions 

The comparison of the measured probe descent spin with the 
calculated spin assessment, as well as the excellent acquisition of 
the communications link to the Orbiter, show that the aerodynamic 
spin design of the probe was entirely adequate. The spin was slow 
enough to provide favorable conditions for timely establishment of 
a communications link with the overflying Orbiter and fast enough 
to provide the science payload sensors a continuous 360-deg survey 
of the local horizon throughout the Jovian descent mission. 
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Nomenclature 

C p = pressure coefficient 

D = maximum payload shroud diameter, m 

d = booster diameter, m 

e = specific energy, J/kg 

F = vector of x -directed fluxes 

G = vector of r-directed fluxes 

H = source vector 

M = Mach number 

p = static pressure, N/m 2 
q = heat flux, W/m 2 

Re = Reynolds number 

t = time, s 

V = conservative variables in vector form 
u, v = velocity components, m/s 
x, r = coordinate directions, m 
p = density, kg/m 3 

cr = stress vector, N/m 2 
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Introduction 

BULBOUS payload shroud is generally employed to accom¬ 
modate an increase in payload volume of the satellite in a 
launch vehicle. A complex flowfield is observed around such a con¬ 
figuration due to the formation of a terminal shock, its interaction 
with the turbulent boundary layer, and flow separation in the boat-tail 
region during ascent flight in the transonic-speed range. The bound¬ 
ary layer in the boat-tail region is considerably thicker; therefore, 
the evaluation of viscous-inviscid interaction becomes necessary. 
Assessment of the separation region and knowledge of the surface 
pressure distribution are required for acoustic insulation design pur¬ 
poses. 

A literature survey shows that experimental and numerical ef¬ 
forts have been made to analyze the complex fluid mechanics of 
the bulbous payload shroud of a launch vehicle. Experimental in¬ 
vestigations were limited to the measurement of surface pressure 
distribution and oil flow patterns. 1 ' 2 Numerical simulations have 
been performed to examine the shock/turbulent boundary-layer in¬ 
teraction and separated flow region in the vicinity of the boat-tail 
region. 3 4 For the various heat shield geometries, the available data 
are not sufficient to determine the movement of the terminal shock, 
the shock-induced separated flow on the forebody cylinder, and the 
separated flow caused by the boat tail. 

The present work employs a computational fluid dynamics ap¬ 
proach to analyze such a complex fluid dynamics problem. The time- 
dependent, compressible, turbulent Reynolds-averaged Navier- 
Stokes equations are solved using a finite volume discretization in 
conjunction with a three-stage Runge-Kutta time-stepping scheme. 
The main purpose of the present numerical simulation is to study 
the formation of the terminal shock and its movement over the fore¬ 
body cylinder in the transonic Mach number range. The behavior 
of the separated flow in the boat-tail region is also investigated in 
the Mach number range 0.8 < M < 3.0 and Reynolds number range 
3.314 x 10 7 ^1.682 x 10 7 /m. A comparison is made with the exper¬ 
imental results obtained from the wind tunnel. 

Governing Equations and Numerical Scheme 

The axisymmetric, compressible, turbulent, Reynolds-averaged 
Navier-Stokes equations can be written as 

dU BF 1 B(rG) H 

— _l_ — -j--= — (1) 

Bt 3x r dr r 

where 

U = [p, pu, pv, pe] T 

F = [pu, pu 2 - o xx , puv - a xr , (pe - o xx )u - a xr v + q x ~f 

G = [pv, puv — a xr , pv 1 — a rr , (pe — a rr )v — a rx u + # r ] 

H = [0, 0, cr + , Of 

Reynolds stress and turbulent heat fluxes are modeled by introduc¬ 
ing an isotropic eddy viscosity term. The coefficient of molecular 
viscosity is calculated according to Sutherland’s law. The closure 
of the system of equations is achieved by introducing the Baldwin- 
Lomax turbulence model. 5 The temperature is related to pressure 
and density by the perfect-gas equation of state. 

The numerical algorithm employs a finite volume discretization 
technique. The computational domain is divided into a number of 
quadrilateral cells. The integral conservation equations are applied 
to each computational cell in turn. The convective fluxes are calcu¬ 
lated at the center of the cell, whereas the viscous fluxes are defined 
at the vertices of the cell. 6 A system of ordinary differential equa¬ 
tions in time is obtained after summing the flux vectors on each face 
of the cell. Temporal integration is performed using the multistage 
time-stepping method. 7 Further details of the numerical procedure 
are given in Ref. 8. 

Four different types of boundary conditions are required for the 
computation of the flowfield, i.e., inflow, solid wall, outflow, and 
centerline conditions. At the inflow, all of the flow variables are 
prescribed at freestream conditions. On the heat shield wall, the no¬ 



slip and isothermal wall conditions are considered. For the subsonic 
flow, nonreflecting far-field boundary conditions are applied at the 
outer boundary of the computational domain. For supersonic flow, 
all of the flow variables are extrapolated at the outflow from the 
vector of conserved variables V. At the line of symmetry, ahead of 
the heat shield, an image cell is imposed on the solved variables. 

The aforementioned numerical procedure is employed to obtain 
the flowfield along a typical bulbous payload shroud. The maximum 
payload shroud diameter of the model is 0.04 m, whereas the booster 
diameter is 0.035 m. The boat-tail angle is 15 deg measured clock¬ 
wise from the axis with reference to the oncoming flow direction. 
For the blunt-nosed cone, the inclination at the forebody junction is 
20 deg and the total length of the shroud from the stagnation point 
to the boat tail is 0.083 m. These subscale dimensions were chosen 
after considering the blockage and compatibility conditions with the 
model system support of the wind tunnel. 1 ' 9 Grid-independent tests 
and the convergence criterion are discussed in Ref. 8. The compu¬ 
tation is performed using 132 x 62 grid points. Figure 1 shows an 
enlarged view of the grid. 

Results and Discussion 

The current experimental data consist of schlieren photo¬ 
graphs and surface pressure distributions. 9,10 The flowfield and re¬ 
sulting wall pressure distributions are computed for the freestream 
Mach number range 0.8 < M < 3.0 and Reynolds number range 
3.314 x 10 7 -4.682 x 10 7 /m based on the wind-tunnel test condi¬ 
tions. 9 Figure 2 shows the comparison between the density con¬ 
tour plots and schlieren photographs for M = 0.8 and 0.9. It can 
be observed from Fig. 2 that all of the essential flowfield features 
of the transonic flow, such as supersonic pocket, terminal shock, 
and expansion and compression regions, are very well captured and 
compare well with the schlieren photographs. The terminal shock 
just reaches the end of the forebody cylinder at M — 0.95, as shown 
in Ref. 8. The density contour plots reveal that the supersonic region 
increases with increasing freestream Mach number, and as a result, 
the terminal shock moves downstream with increasing freestream 
Mach number. Comparison between the shock locations of the nu¬ 
merical and experimental results is shown in Fig. 3. It can be seen that 
the terminal shock crosses the forebody cylinder at about M = 0.97. 
The location of the terminal shock is found to be a nonlinear function 
of Mach number. 

In recent numerical studies, 8 a local flow separation is observed in 
the boat-tail region and is found to be a function of Mach number. A 
schematic sketch of the separated flow in the boat-tail region at the 
transonic and supersonic Mach numbers, which are reconstructed 
from the velocity vector 8 and density contour plots, is shown in 
Fig. 4. For the transonic case, a shear layer is formed, and it con¬ 
tains a recirculating region as shown in Fig. 4a. The downstream 
boundary layer is found to be thick, which is nearly the boat-tail 
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Fig. 2 Comparison between density contours and schlieren picture. 



X/D 

Fig. 3 Shock location at various freestream Mach numbers. 

height. Figure 4b shows flowfield features for the supersonic case. 
The flow turns downstream of the boat tail through an expansion 
fan with a resulting mixing layer that separates the supersonic fluid 
above the subsonic recirculating region below it. The pressure in 
the recirculating region is generally uniform and below that of the 
freestream pressure. As the mixing layer impinges on the down¬ 
stream wall, the flow turns back toward the freestream direction and 
a recompression shock wave is typically formed. 

The surface pressure coefficient C p vs nondimensional length 
x/D of the heat shield for various values of freestream Mach num¬ 
bers is presented in Figs. 5 and 6. In the case of the transonic-speed 
range, the computed results compare well with the experimental re¬ 
sults along the entire length of the shroud except at the cone-cylinder 
junction in the transonic range, where the maximum expansion oc¬ 
curs. Near the stagnation point and on the spherical region, a steep 
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Fig. 4 High-speed flow over boat-tail region. 



Fig. 5 Pressure distribution along heat shield; M < 1.0. 


fall in the pressure coefficient was noticed, which is partially recov¬ 
ered along the cone. The experimental values also agree with this 
trend from the station x/D = 0.22 onward, where the first pressure 
measurement was made. On the forebody cylinder, the pressure dis¬ 
tribution pattern also confirms the movement of the terminal shock, 
as observed in the density contour plots. In the boat-tail region, the 
expansions and compressions are comparatively mild. In the case 
of the supersonic Mach number range, the expansion at the cone 
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Fig. 6 Pressure distribution along heat shield; M < 3.0. 

forebody and the boat-tail region reduces with increasing super¬ 
sonic Mach number. 

Conclusions 

A numerical experiment of axisymmetric turbulent viscous flow 
over a bulbous heat shield is performed by employing a three-stage 
Runge-Kutta time-stepping scheme. Turbulence closure is achieved 
using the Baldwin-Lomax turbulence model. The flowfield visual¬ 
ization of the terminal shock and the separated region helps in a 
systematic understanding of flow structure under various freestream 
Mach numbers. The following observations are made based on the 
described numerical simulations. 

1) The terminal shock moves downstream with increasing 
freestream Mach number. The strength of the terminal shock initially 
increases with Mach number and later decreases. The location of 
the terminal shock is found as a nonlinear function of the freestream 
Mach number. 

2) The separation zone in the boat-tail region is found as a function 
of the freestream Mach number. The flowfield features of the sepa¬ 
rated region are found to be different in the transonic and supersonic 
Mach number ranges. 
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Nomenclature 

a = speed of sound, m/s 
D = nozzle external diameter, m 

F = convective flux vector in x direction, nondimensional 

G = convective flux vector in y direction, nondimensional 

H* = axisymmetric source vector, nondimensional 

K = node number, nondimensional 

M = Mach number, nondimensional 

P = pressure, nondimensionalized with P x 

R = residual, nondimensional/radius, m 

r = radial distance, m 

S = control volume surface, nondimensionalized with D 2 
t = time, nondimensionalized with a^ and D 
U = solution vector, nondimensional 
V = velocity vector, nondimensionalized with a x 
x,y = two-dimensional or axially symmetric coordinates, 
nondimensionalized with D 
y = ratio of specific heats, nondimensional 
S2 = control volume, nondimensionalized with D 3 

Subscripts 

j = jet 

oo = infinity 

Introduction 

N exhaust plume expansion into a supersonic external free¬ 
stream is a complex phenomenon. The plume boundary be¬ 
haves as a curved rigid body, and it causes a curved shock wave. 
Very large gradients exist around these regions. 1 Some experimental 
studies have been performed to obtain plume behavior at different 
jet exit pressure ratios. 2-4 There are several methods for predicting 
exhaust plume boundaries. 12 ' 5 ' 6 

In this study, axisymmetric Euler equations were solved by using 
a finite volume technique on unstructured triangular grids. Jet plume 
boundaries and other flowfield details were obtained for supersonic 
axisymmetric jets in a parallel supersonic external flow. The jet-to- 
freestream pressure ratio varied between 1 and 1500. For very high 
pressure ratios, computations are started with a low pressure ratio, 
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